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We propose a method to produce entangled states of several par- 
ticles starting from a Bose-Einstein condensate. In the proposal, a 
single fast n/2 pulse is applied to the atoms and due to the collisional 
interaction, the subsequent free time evolution creates an entangled 
state involving all atoms in the condensate. The created entangled 
state is a spin-squeezed state which could be used to improve the 
sensitivity of atomic clocks. 

The possibility of creating and manipulating entangled 
states of many-particle systems has recently boosted the field 
of quantum information since it may yield new applications 
which rely on the basic principles of Quantum Mechanics!! 
So far, up to four atoms have been entangled in a controlled 
wayBEl. The. experimental achievement of Bose-Einstein 
condensationtaEl has also raised a lot of attention since it 
may lead to applications in several fields of Science. Some 
of these applications are based on the fact that condensates 
can be considered as pure states at the single particle level, 
which is a crucial requirement for the production of entangled 
state. Thus, a natural question is to investigate whether Bose- 
Einstein condensates can also be used in some applications of 
quantum information. In this paper we show that this is indeed 
the case. We propose a method to obtain substantial entangle- 
ment of a large number of atoms with present technology. In 
our proposal, a single resonant pulse is applied to all atoms 
in the condensate, and the collisional interaction entangles the 
atoms in the subsequent free evolution. 

Consider a set of N two-level atoms confined by some 
external trap. In order to describe the internal properties of 
these atoms, it is convenient to let the internal states \a) n 
and \b) n of the n-th atom represent the two states of a fic- 
titious spin 1/2 particle with angular momentum operators 
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l/2(\a)(a\ n -\b)(b\ n ),ji n) = l/2(\b)(a\ n + \a)(b\ 



and jy — i/2(\b)(a\ n — \a)(b\ n ). We consider collective 
effects of the atoms which are described by total angular mo- 
mentum operators, J = Yln=i3^ '• The entanglement prop- 
erties of the atoms can be expressed in terms of the variances 
and expectation values of these operators. In the appendix we 
show that if 



a = N(AJ Hl ) 2 



< 1, 



(1) 



where = n ■ J and the n's are mutually orthogonal unit 
vectors, then the state of the atoms is non-separable (i.e. en- 
tangled). The parameter £ 2 thus characterizes the atomic en- 
tanglement, and states_with £ 2 < 1 are often referred to as 
"spin squeezed states'13. In the following we show how to re- 
duce £ 2 by several orders of magnitude using the collisional 
interactions between atoms in a Bose-Einstein condensate. 
We consider a two component weakly interacting 



Bose-Einstein, condensate, which is present in several 
laboratoriesQu, and we assume that the interactions do not 
change the internal state of the atoms. This situation is de- 
scribed by the second quantized Hamiltonian 



H = 



j=a,b 
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U ab I d 3 rtft(r)*t(r)* a (r)tf 6 (r), 
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where H$ j is the one particle Hamiltonian for atoms in state 
j including the kinetic energy and the external trapping po- 
tential Vj(r), $!j(r) is the field operator for atoms in the state 
j, Ujk — At:h 2 ajk/ra is the strength of the interaction be- 
tween particles of type j and k, parameterized by the scatter- 
ing length djk, and m is the atomic mass. 

Assume that we start with a Bose-Einstein condensate in 
state | a) at very low temperature (T ~ 0), so that all the 
atoms are in a single particle (motional) state \<fro). A fast 
7r/2 pulse between the states \a) and \b) prepares the atoms in 
the state l^o)®^ <8> (|a) + \b))® N /2 N ' 2 which is an eigenstate 
of the J x operator with eigenvalue N/2. If we choose in ([j]) 
n\ = [0, cos(0), sin(0)] and ni along the x axis of the ficti- 
tious angular momentum, we have Q = 1 at t = 0. By using 
the equations of motion of the angular momentum operators 
in the Heisenberg picture we find the time derivative of Q at 
t = Q 
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(N-l)(U aa + U bb -2U ab ) 
2h 



d 3 r\ 
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This equation immediately shows that spin squeezing will be 
produced for certain, angles 9 if U aa + U bb ^ 2U ab as it is in 
the Na experimentsH. 

To quantify the amount of squeezing which may be ob- 
tained, we assume identical trapping potentials V a (r) = 
Vb(r) and identical coupling constants for interaction between 
atoms in the same internal state a aa = a bb > a ab . Physi- 
cally, this could correspond to the \F — l,Mp — ±1) hy- 
perfine states of Na trapped in an optical dipole trap. Due 
to the symmetry of these states their scattering lengths and 
trapping potentials will be identical; moreover, due to angular 
momentum conservation there are no spin exchanging colli- 
sions between these states as required by our model. This is 
exactly the experimental setup used in Ref. ^[ In this exper- 
iment it is shown that these states have an antiferromagnetic 
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interaction a a b < a aa which according to enables the pro- 
duction of squeezed states. To avoid spin changing collisions 
that populate the state \Mp = 0)EJ one has to slightly mod- 
ify such an experimental set-up. If the F = 1 manifold is 
coupled to the F = 2 with a far off-resonant blue detuned 
7r-polarized microwave field, the \F — l,Mp = 0) state is 
raised in energy with respect to the \F = 1, Mf = ±1) states, 
since the Clebsch-Gordan coefficient for the \F = 1, Mp = 
0) — > \F = 2,Mp = 0} is larger than the coefficients for 
the \F = 1, Mp = ±1} -> \F = 2, Mp = ±1) transitions, 
and spin exchanging collisions become energetically forbid- 
den. If for instance one chooses a detuning 6 = (27r)25 MHz 
and a resonant Rabi frequency for the 1 — * 1 transition of 
Q = (27r)2 MHz, the energy separation is AE = 640nK. 
With a typical chemical potential fx « 220 nKEJ, this energy 
separation is much higher than the available energy in the col- 
lisions and the \F = 1, M = 0) state is completely decoupled. 



culated exactl yB In the limit of large N, the minimum obtain- 
able squeezing parameter is = \(jj) 2 ^ , which indicates 
that our proposal might produce a reduction of £ 2 by a factor 
of ~ TV 2 / 3 which would be more than three orders of magni- 
tude with 10 5 atoms in the condensate. 

In contrast to the simplified Hamiltonian H spin the real 
Hamiltonian (|^) will also entangle the internal and motional 
states of the atoms which is a source of decoherence for the 
spin squeezing. To quantify this effect we have performed 
a direct numerical integration following the procedure devel- 
oped in Ref. O, We split the whole Hilbert space into orthog- 
onal subspaces containing a fixed number of particles N a and 
Nb = N — N a in each of the internal states, respectively. In 
each subspace we make a Hartree-Fock variational ansatz in 
terms of three-dimensional spatial wavefunctions (f> a (N a ,t) 
and 0h(iVb,t), which are evolved according to the time de- 
pendent coupled Gross-Pitaevskii equations. This is an ap- 
proximation to the full problem which is valid in the limit 
of low temperatures and short times, where the population of 
the Bogoliubov modes is small. Particularly it is a good ap- 
proximation in our case with a a b < a aa , abb where there are 
no demixing instabilities. With this procedure, the decoher- 
ence induced by the entanglement with the motional state is 
effectively taken into account. Together with the prediction 
from the simple Hamiltonian H spin the result of the simula- 
tion is shown in Fig. [j]. The two curves are roughly in agree- 
ment confirming that the system is able to approximate the re- 
sults of the Hamiltonian H sp i n . The numerical solution shows 
fluctuations due to the oscillations of the spatial wavefunc- 
tion. The large dips at tot « 4, 9, 13, and 18 are the points 
where the atomic cloud reach the initial width. At these in- 
stants the overlap of the wavefunctions is maximal and the 
two curves are very close (up to a factor of two). The small 
dips at u)t =2, 7, 11, and 16 corresponds to the points of max- 
imum compression. With the realistic parameters used in the 
figure, our simulation suggests that three orders of magnitude 
squeezing is possible. Also, note the time scale in the fig- 
ure. The maximally squeezed state is reached after approx- 
imately two oscillation periods in the trap. For a fixed ratio 
of the scattering lengths a a b/a aa , the optimal time scales as 

(a aa /do)~ 2/5 ^V~ 1/15 , where d = y/KJJmtJ) is the width of 
the ground state of the harmonic potential. 



In Ref. 10 a much smaller energy difference is shown to ex- 
clude spin exchanging collision and we therefore expect that 
much weaker fields will suffice. 

The assumption a aa = abb has several advantages. Firstly, 
it reduces the effect of fluctuations in the total particle number. 
If CLaa 7^ o-bb the mean spin performs a N dependent rotation 
around the z-axis, and fluctuations in the number of particles 
introduces an uncertainty in the direction of the spin which ef- 
fectively reduces the average value and introduces noise into 
the system. With a aa = a^ the mean spin remains in the 
x-direction independent of the number of atoms in the trap. 
Secondly, this condition ensures a large spatial overlap of dif- 
ferent components of the wavefunction. After the ir/2 pulse 
the spatial wavefunction is no longer in the equilibrium state. 
That is, due to the atomic repulsions (which are now different 
than before since the atoms are in different internal states), the 
spatial distribution of the atomic cloud will start oscillating. 
Furthermore, since the state of the system is now distributed 
over a range of number of particles in the \a) state (N a ), and 
since this number enters into the time evolution, the N a de- 
pendent wavefunctions <fi a and <j>b are different for particles in 
the states \a) and \b). With a aa — abb, 4>a and <j>b are identical 
if N a equals the average number N/2. In the limit of large N, 
the width of the distribution on different N a 's is much smaller 
than N a and all the spatial wavefunctions are approximately 
identical <fi a (N a , t) « <j>b{N a ,t) w (f>o(t). This relation is 
only true if a aa > a a b where small deviations from the aver- 
age wavefunction perform small oscillations. In the opposite 
case the deviations grow exponentiallyU resulting in a reduc- 
tion of the overlap of the a and the b wavefunctions and a 
reduced squeezing. The advantages mentioned above could 
also be achieved with a aa ^ cibb by using the breath-together 
solutions proposed in Ref. [lX|. 

Before analyzing quantitatively the complete system, we 
estimate the amount of spin squeezing we can reach with our 
proposal by using a simple model. Assuming the same wave- 
function <j>o for both | a) and |6) atoms is constant and inde- 
pendent of the number N a , the spin dependent part of the 
Hamiltonian (||) may be written as H spin = h\J^, where \ 
depends on the scattering lengths and the wavefunction tfio. 
The spin squeezing produced by this Hamiltonian can be cal- 
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FIG. 1. Reduction in the squeezing parameter £ 2 . A fast it/ 2 
pulse between two internal states is applied to all atoms in the con- 
densate. The subsequent free time evolution results in a strong 
squeezing of the total spin. The angle 9 is chosen such that £f is 
minimal. The solid line is the results of a numerical integration 
(see text). For numerical convenience we have assumed a spheri- 
cally symmetric potential V(r) — rnuj 2 r 2 /2. The parameters are 
a-aa/do = 6 * 1CT 3 , a bb = 2a ab = a aa , and N = 10 5 . The 
dashed curve shows the squeezing obtained from the Hamiltonian 
H Bp i n — h\Jz- The parameter % is chosen such that the reduction 
of J x obtained from the solution in Ref. ^ is consistent with the re- 
sults of Ref. |ll[ 

The analysis so far has left out a number of possible im- 
perfections. Specifically, we have assumed that all scattering 
length are real so that no atoms are lost from the trap and we 
have not considered the role of thermal particles. To estimate 
the effect particle losses we have performed a Monte Carlo 
simulatior£3 of the evolution of squeezing from the Hamil- 
tonian H sp i n . The particle loss is phenomenologically taken 
into account by introducing a loss rate T which is identical 
for atoms in the |o) and the \b) state. In Fig. ^ we show the 
obtainable squeezing in the presence of loss. Approximately 
10 % of the atoms are lost at the time %t « 6 * 1CP 4 where 
the squeezing is maximally without loss. With the parame- 
ters of Fig. [y this time corresponds to roughly two trapping 
periods. Such a large loss is an exaggeration of the loss com- 
pared to current experiments and the simulation indicate that 
even under these conditions, squeezing of nearly two orders of 
magnitude may be obtained. On the other hand, the effects of 
thermal particles can be suppressed at sufficiently low temper- 
atures but due to the robustness with respect to particle loses 
shown in Fig. 2, we expect to obtain high squeezing even at 
some finite temperatures. 

In conclusion we believe that we have presented a sim- 
ple and robust method to produce entangled states of a 
large number of atoms with present technology. The pro- 
duced entangled state are interesting in fundamental physics 
and they also fraye possible technological applications in 
atomic clocksli3o, where the projection noise (A Jg) 2 is cur- 
rently the-main source of noiset3. In particular Boyer and 
Kasevichli-3 have shown how to use an entangled state with 
half of the particles N/2 in one internal state a, and the other 



half in the other state b to obtain the Heisenberg limit in atomic 
interferometric measurements. On the other hand, it has been 
shown that if the atoms are prepared in a state with £ 2 < 1 
before they are injected into an atomic clock, one can reduce 
the frequency noise (variance in the frequency measurements) 
or the measuring time to obtain a desired precision by a factor 
£ 2 as comparedip the case in which one uses atoms in an un- 
corrected stateEE If the squeezing produced by our proposal 
is transferred into a suitable clock transition, and the atomic 
cloud is allowed to expand so that the role of collisions is re- 
duced during the frequency measurement, the states could be 
directly applied in the current set-up of atomic clocks. Other 
theoretical proposals few noise reduction in neutral atom sys- 
tems have been madeE3, and a weak squeezing, of the spin 
has recently been produced experimentallyLZfEJ. However, 
our proposal has the considerable advantage that it offers a 
very strong noise reduction, and it is directly applicable to ex- 
isting experimental set-ups. In future experiments with neg- 
ligible particle loss even for very long interaction times, the 
Hamiltonian H sp i n could also be used tanroduce maximally 
entangled state of any number of atomsElo, which could re- 
duce the frpauency noise to the fundamental limit of quantum 
mechanics! 2 ^. 
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FIG. 2. Quantum Monte Carlo simulation of squeezing in the 
presence of loss. The full line shows squeezing obtained from the 
Hamiltonian H sp i n with particle loss described by a constant loss 
rate V = 200x- The dashed curve shows squeezing without particle 
loss. 

Appendix 

Here we present the derivation of Eq. ([[]) as a criterion for 
entanglement. An iV-particle density matrix p is defined to be 
separable (non-entangled) if it can be decomposed into 

k 

where the coefficients are positive real numbers fulfilling 
E/c Pk = 1, and pi' is a density matrix for the i'th parti- 
cle. The variance of J z may be described as (A J z ) 2 = -j — 

E fe Pk Ei(#)fc + E* PMl - (J*?, and using Cauchy- 
Schwarz's inequality and + (j^) 2 + (ji l) ) 2 < 1/4 
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we find three inequalities for separable states ^2 k Pk{J z )k — 

(■h)\ -E k PkEi(j^)l > -f + E k PkEM l) )l + 

and (J x ) 2 < NY, k Pk^{3^)l From these in- 
equalities we immediately find that any separable state obeys 
£ 2 > 1 and hence any state with £ 2 < 1 is non-separable. 
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